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Abstract. The next generation of gravitational wave detectors holds out the prospect 
of detecting a stochastic gravitational wave background generated in the very early 
universe. In this article, we review the various cosmological processes which can lead to 
^ ^ such a background, including quantum fluctuations during inflation, bubble collisions 

\^ ' in a first-order phase transition and the decay of a network of cosmic strings. We 

0^ ■ conclude that signals from strongly first-order phase transitions, possibly at the end 

of inflation, and networks of local cosmic strings are within the sensitivity of proposed 
detectors. However, backgrounds from standard slow-roll inflation and the electroweak 

■ phase transition are too weak. 

Pi 1. Introduction 

■ An outstanding achievement of modern cosmology has been the detection of anisotropies 
in the cosmic microwave background These anisotropies provide a snapshot 

^ ' of the universe about 400,000 years after the Hot Big Bang, just as the universe 

■ became transparent to electromagnetic radiation. Despite this observational triumph, 
as yet it has been insufficient to differentiate between competing paradigms for galaxy 
formation, nor has it shed much light on cosmological processes taking place before 
photon decoupling. There are, however, other types of radiation, such as gravitational 
radiation, which penetrate through this electromagnetic surface of last scattering, 
travelling virtually unaffected since their emission, even when this occurred in the first 
few fractions of a second. Of course, for gravitons this remarkable transparency is due to 
their very weak interactions with ordinary matter which, in turn, makes them difficult 
to observe. However, pioneering experiments have been proposed which could detect 
stochastic backgrounds of gravitational waves generated in the early universe over a 
range of frequencies. 

Gravitational radiation detected with a particular frequency today would have been 
created at a characteristic time (if generated by classical, causal mechanisms). Assuming 
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that the radiation is emitted at a time te before equal matter-radiation (tg < teq ~ 40, 000 
years) and that it is created with a wavelength comparable to the horizon A(tc) ~ te, then 
its frequency today is given by / ~ z^^{tcqtc)~^^'^ where the redshift Zeq ~ 2.3 x lO'^Vloh'^ 
and 'little /i' is the rescaled Hubble constant lying in the range 0.4 < h < 0.9. The 
quantity that is measured by a particular experiment is the dimensionless amplitude 



where fig{f) is the gravitational radiation density contribution to the universe's density 
per octave at a frequency /, and pc is the critical energy density. Clearly, from (|l|) 
gravitational radiation with higher frequencies will be more difficult to detect. 

There are currently four frequency bands available for studying gravitational 
radiation 0. First, there are the microwave background anisotropies created by tensor 
modes which may have been detected by the COBE-DMR experiment. In inflationary 
scenarios, the precise signal ratio due to tensor versus scalar modes is sensitively model- 
dependent, but it could be as high as 50%. Assuming that the entire COBE signal 
is due to gravitational waves implies a weak upper bound of f2g < 7 x 10^^^ at a 
frequency of 3 x 10~^^/iHz. A second upper limit on f2g comes from pulsars, since 
a stochastic background would lead to timing noise in measurements of the incoming 
periodic electromagnetic signal. A recent limit is Qgh'^ < 6 x 10~^ at 4 x 10~^Hz ||^, 
though the statistical veracity of this result has been questioned with re-analyses of 
the same data suggesting both stronger and weaker bounds 0, H. Thirdly, proposed 
ground based detectors will study frequencies around lOOHz with a maximum sensitivity 
of about Qgh"^ ^ 10"'^ for the first generation of LIGO detectors and Qgh"^ ^ 10~^° for 
the next generation LIGO and VIRGO detectors. Finally, the proposed space-based 
interferometer LISA will have a sensitivity of Qgh^ ^ 10~^° at about 10~^Hz. 

2. Slow-roll inflation 

The inflationary paradigm ||^, ^ ^ is strongly motivated because it resolves a number 
of shortcomings of the standard Hot Big Bang cosmology, including the horizon, 
flatness and monopole problems. However, possibly the most significant testable 
prediction of inflationary scenarios is a nearly scale-invariant spectrum of adiabatic 
density perturbations produced through quantum mechanical effects when the very 
early universe is dominated by the potential energy of some scalar inflaton field (p. 
An important by-product of this process is the excitation of tensor modes resulting in 
a background of gravitational waves. The spectra for scalar and tensor modes created 
in a slow-roll inflationary model are the following. 
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where H{(f)) is the Hubble parameter during inflation, 0^ is the value of the inflaton field 
when the mode with comoving wavenumber uj = 2Txf leaves the Hubble radius H~^, and 
the slow-roll parameter is 

with the prime ' denoting differentiation with respect 0. Here, e(0) is related to the 
speed at which rolls down the potential ^(0) (less than unity during inflation). Once 
a particular mode is created, it is driven outside the Hubble radius (or 'horizon') by the 
rapid expansion and it effectively 'freezes' until it returns inside the horizon during the 
subsequent radiation and matter dominated eras. Lengthscales corresponding to the 
size of the observed universe went outside the horizon early when the inflaton field had 
a value denoted by 06o; after this there were about 60 e- foldings of expansion before the 
end of inflation when e(0end) = 1- 

For a specific inflaton potential 1^(0), the dynamics of the Hubble parameter in the 
Hamilton- Jacobi formalism is given by |^ 

^,,_l^^,_32^ (4) 

m^i m^i 47r 

The Hubble parameter monotonically decreases during slow-roll inflation, so from (|^) the 
largest contribution to the gravitational wave spectrum is due to modes that were driven 
outside the horizon early in inflation and have just come back inside the horizon at the 
present day. These scales correspond to those probed by the COBE-DMR experiment 
and the observed anisotropies can be used with (|^) to normalize Hqq = H^cpQo) . Given 
this initial condition, (^ can be used to calculate the Hubble parameter during inflation 
and hence the resulting spectrum of gravitational waves. Paradoxically, as we shall 
see, the larger the contribution of a model to the COBE signal, the smaller the signal 
will be in the frequency bands corresponding to LISA and LIGO/VIRGO; the Hubble 
parameter must decrease more rapidly in order to maintain the condition that there be 
60 e- foldings before the end of inflation |TD[. 

In order to illustrate this point, we shall calculate the spectrum of gravitational 
radiation produced by an inflationary model similar to polynomial chaotic inflation, but 
slightly modified to make analytic calculations tractable Using a simple method 
to obtain the COBE normalization the potential for this model is given by 



.,.)^a.rx.o-.X,(:.^)-'(|-y(:-£4|, ,) 



where 06o = — (a"^pi/4A/7r)(l + 240/q;)^/^ and we require a > 1. For this model, the 
slow-roll and Hubble parameters are given by 

167r0^ mpi V a J \06o/ 
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For oj > Ueq, substituting into (|2D and appropriately redshifting yields 



a J 



(7) 



/ 24()\ r 4 / 

fi.(cu)/i2 ^ 4.6 X lO^^M 1 + — 1 log 

^ V a / [ 240 + a \aoHo 

For uj < Ueq, that is, scales which come inside the horizon later during the matter era, 

we have less redshifting and (0) must be multiplied by an additional factor (cUeq/cu)^. 

To delineate the predictions of different slow-roll inflation models, one can use 
different values of a; the obvious limiting cases being a = 2, the simplest chaotic model, 
and a = 80 (say) for which almost all the COBE signal is due to gravitational waves 
(this case corresponds closely to power law inflation with an exponential potential). The 
results of these cases are plotted in figure 1. The a = 2 model has a low contribution 
to COBE but the spectrum for u > c^eq decreases slowly down to an amplitude of 
flgh"^ ~ 10^^^ at a few Hz. In contrast, the model with a = 80 has a much larger 
contribution to COBE, but has a very small amplitude for scales uj > uJcq] in the LISA 
frequency band this has already fallen to Qgh"^ ^ 10"^^, eight orders of magnitude 
smaller than for a = 2. Nonetheless, both these models produce signals well below the 
sensitivity of proposed interferometers by several orders of magnitude. 

The largest possible COBE normalization provides a weak upper limit on the 
amplitude of gravitational waves for all scales u > ujeq which must be below the level 
Qgh"^ ^ 10~^^. It may be possible to carefully construct a slow-roll potential which 
would almost achieve this limit in the LISA and LIGO/VIRGO frequency bands, but 
undoubtedly the potential would look rather contrived. In any case, it is clear that 
the prospects for a detectable slow-roll signal are bleak, so we shall have to consider 
less orthodox inflationary models. In the next section, we consider extended and hybrid 
inflation, but there are also more speculative superstring-inspired models which can give 
a larger gravitational wave background (for example see |1l2| ). 



3. First-order phase transitions 

Cosmological phase transitions in which symmetries are spontaneously broken are now 
an integral part of modern cosmology. In a first-order phase transition the field 
becomes trapped in a metastable local minimum of the potential — the false vacuum. 
The transition to the true vacuum takes place by the nucleation of bubbles and their 
subsequent rapid growth. When these vacuum bubbles collide copious amounts of 
gravitational radiation can be emitted, particularly if the relative velocity of the walls 
is relativistic as in strongly first-order phase transitions. It has been estimated that 
the maximum contribution to the gravitational wave background from such a transition 
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Figure 1. Summary of the potential cosmological sources of a stochastic gravitational 
radiation background, including inflationary models, first-order phase transitions and 
cosmic strings, as well as a primordial 0.9K black-body graviton spectrum (the 
analogue of the black-body photon radiation). Also plotted are the relevant constraints 
from the COBE measurements, pulsar timings, and the sensitivities of the proposed 
interferometers. Notice that local cosmic strings and strongly first-order phase 
transitions may produce detectable backgrounds, in contrast to standard slow-roll 
inflation models. 

at a frequency of 

/„„«3X10-H.(|)(^)""(^), (9) 

where F = Fq exp(/5t) is the nucleation rate of bubbles, is the temperature of the phase 
transition, is the relevant Hubble parameter and jV* is the number of relativistic 
degrees of freedom. Assuming that the electroweak phase transition is strongly first- 
order imphes Qgh'^ 10~^ at a frequency of /max ~ 10~^, inside the LISA band with a 
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detectable amplitude, as illustrated in figure 1. However, the minimal standard model is 
currently believed to have only a weakly first-order transition, if not second-order, with 
the bubble walls reaching velocities well below the speed of light. In this case, Qgh"^ 
has been estimated to be around 10~^^ which is well outside LISA's sensitivity M . 
Nevertheless, there are a number of well-motivated extensions to the standard model 
(like supersymmetry) which could entail symmetry-breaking just above the electroweak 
scale; if strongly first-order, such phase transitions would create a distinctive and 
detectable LISA signal. Fortuitously from this point of view, LISA has a very interesting 
frequency response range. 

While the prospects for detecting gravitational waves from slow-roll inflation seem 
poor, there are a number of inflationary models which do not end in the standard re- 



heating scenario. These include extended inflation [|15l and hybrid inflation [|T6[ which 
exit through a phase transition; this provides an extra, potentially more powerful source 
of gravitational waves if this flnal phase transition is strongly flrst-order [|l7l. The 
radiation could be detected by the advanced LIGO and VIRGO detectors if the re-heat 
temperature of the universe was close to lO^GeV or lO^GeV — scales which are far from 
ruled out by current observations. We have included a sample spectra for such a model 
in flgure 1. 



4. Cosmic string networks 

Cosmic strings are line-like topological defects which may form during phase transitions 
in the early universe (see, for example, [|l8|). Due to the large string mass per unit 
length (typically about lO^^gcm"^ for GUT scale strings), they may have acted as the 
initial seeds for the formation of large-scale structure [|l^, In general, a network of 
strings will evolve towards a self-similar scaling solution by the production of loops and 
the subsequent emission of radiation into a preferred channel. For local strings, this 
channel is gravitational radiation, thus creating a stochastic background of radiation 
over frequencies from lO^^Hz to 10~^^Hz. 

Scale-invariant string evolution implies that the density of the string network 
Poo = fJ'C/t'^ remains constant relative to the background density. Statistically, the 
network has the same properties at any two times, except for a universal scaling with 
respect to the growing horizon size, that is, there exists a characteristic length scale 
L = C,~^/'^t. This scaling regime has been shown to exist numerically [^, ^ and 



various parameters have been estimated, notably C ~ 13 during the radiation era. 

The radiative dynamics of the strings is a crucial factor determining the precise 
contributions at various frequencies. The power radiated by a string loop is normally 
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written as a sum over loop harmonics that is, 

P = Y.Pn = TG/^' (10) 



n=l 



where (F) ^ 65 |^ is a constant parametrized solely by the particular loop trajectory 
and not its length, /i is the string mass per unit length and n^, is some cut-off introduced 



by the effects of radiation backreaction pSf (see also p6 



Assuming that M is constant and for reasonable values of ra^,, one can deduce that 
the spectrum of gravitational radiation due to loops created during the radiation era 
is independent of frequency and only very weakly dependent on the loop radiation 
spectrum. The amplitude of this flat spectrum is given by 

where F{x) = ((1 + x)^^"^ — l)/x, fij. = Pr/pc, Pr = 3/327rGt^ and a is the constant loop 
production size with respect to the horizon. In contrast, the radiation background 
today due to loops produced in the matter era is highly dependent on the loop 
radiation spectrum. Assuming that P„ oc one can deduce that ^g{uj) oc uj^~'^ —>■ 
logi^lfl g{uj)] = A + {1 — q) logio[c^]- Therefore, the amount of radiation produced in the 
matter era which feeds into higher frequencies is extremely sensitive to q and n*. It is 
precisely these frequencies that are relevant for the pulsar timing experiment p7 . 



One can investigate these effects more fully as in ref. using a numerical algorithm 
which calculates the spectrum of radiation for given spectra and cosmological parameters 
p7| . The parameters that were used are Qq = 1, h = 0.5 and the electroweak model 
particle spectrum with three lepton generations. It was found that the contribution to 
the pulsar timing frequency was the same for all g > 2 irrespective of the value of n*, 
while for g < 2 it was sensitive to the value of n^:. In particular, for g = 4/3 and n* = oo 
the amplitude of the spectrum will be in conflict with the constraint from pulsar timings 
for the string energy density Gp ~ 10~^ normalised to COBE ||2^. However, if g = 4/3 
and < 1000 then the contribution to the pulsar timing frequency will be the same as 
for g = 2 and ra* = oo. The spectrum of radiation emitted by loops is hkely to be cut off 
by the effects of backreaction as shown numerically for Goldstone boson radiation from 
global cosmic strings p6[, therefore it seems sensible to exclude the g = 4/3, large n^, 
spectra on physical grounds (as we discuss in ref. p5|). Using well-motivated values of g 
and n^:, one can use the pulsar timing constraint to deduce that Gp < 5.4(±0.8) x 10~^ 



25| . Using f2o < 1 also reduces the contribution to the pulsar timing frequency from 



matter era loops due to curvature domination at tc ~ ^oto, though we expect this effect 
to be small. 

The contribution from cosmic strings to the LISA and LIGO/VIRGO frequency 
ranges is likely to be slightly lower than that relevant for pulsar timing because particle 
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mass thresholds cause entropy transfers which reduce the relative density contribution 
of any pre-existing decoupled radiation. The cosmic string spectrum shown in figure 
1 illustrates this effect with a gradual rise in the otherwise fiat spectra at frequencies 
around 10~^Hz. This particular rise is caused by particle annihilation near the QCD 
and electroweak phase transitions. If the particle physics model has more degrees of 
freedom at higher energies there could be other steps related to other phase transitions. 
However, the dependence on M is reasonably weak and therefore we can conservatively 
estimate that Vtji^ > 2.0 x lO^^ at IQ-^Hz and n^h"^ > 1.0 x IQ-^ at lOOHz for 
Gn = 1.0 X 10^^. We note also that very precise determinations of the stochastic 
background from cosmic strings at different frequencies would measure A/" in different 
cosmological epochs, providing fascinating insight into the particle content of the early 
universe (at times as early as 10~^^sec for LIGO and VIRGO). 

Finally, we should note that there are other types of cosmic strings which do not 
produce a large background of gravitational radiation, such as those formed when global 
symmetries are broken. The dynamics of global strings is dominated by the production 
of Goldstone bosons and a simple calculation suggests that gravitational waves will be 
suppressed by approximately four orders of magnitude; the global string spectrum is 
also illustrated in figure 1. 



5. Conclusion 



Gravitational waves potentially provide a unique and unparallelled probe of the very 
early universe. The proposed interferometers, both terrestrial and space-based, could 
rule out or severely constrain a number of viable theoretical models, most notably the 
cosmic string scenario. On the other hand, the detection of a primordial background of 
gravitational waves would have a profound impact on our understanding of high energy 
physics and cosmology, providing an unprecedented view of the earliest moments after 
the creation of the Universe. 
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